NACA TN 2QG3 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 


TECHNICAL NOTE 2063 


A COMPARISON OF THEORETICAL AND EXPERIMENTAL WING 
BENDING MOMENTS DURING SEAPLANE LANDINGS 

> 

By Kenneth F. Merten, Jos^ L. Rodriguez, 
and Edgar B. Beck 

Langley Aeronautical Laboratory 
Langley Air Force Base, Viar" 


ROOM 


\NACA^ 


Washington 
April 1950 


MAY 6 1993 


LANGLEY RESEARCH CENTER 
LIBRARY NASA 
HAMPTON, VIRGINIA 


i 




NASA Technical Ubn 



76 01433 4040 


NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 2063 


A COMPARISON OF THEORETICAL AND EXPERIMENTAL WING 

BENDING MOMENTS DURING SEAPLANE LANDINGS 

By Eenneth F. Merten, Jos£ L. Rodriguez, 
and Edgar B. Beet 


SUMMARY 


'A smooth-vat er— landing investigation was conducted with a small 
seaplane to obtain experimental wing— bending— moment time histories 
together with time histories of the various parameters necessary for the 
prediction of wing bending moments during hydrodynamic impact. The 
experimental results were compared with calculated results which include 
inertia— load effects and the effects of air— load variation during 
impact. The responses of the fundamental mode were calculated with the 
use of the measured hydrodynamic forcing functions. From these responses, 
the wing bending moments due to the hydrodynamic load were calculated 
according to the procedure given in R. 4M, No. 2221. This comparison 
of the time histories of the experimental and calculated wing bending 
moments showed good agreement both In phase relationship of the' oscil- 
lations and in numerical values. 

The effects of structural flexibility on the wing bending moment 
were large, the dynamic component of the total moment being as much as 
97 percent of the static component . Changes in the wing bending moment 
due to the variation in air load during impact were of about the same 
magnitude as the static water— load component. 


INTRODUCTION 


Recent trends In the design of aircraft . have led to an important 
increase of the stresses produced in wings by landing impacts. Two 
significant factors contributing to these Increased stresses during 
landing are an Increased proportion of the seaplane weight in the wings 
and an increased structural flexibility, since, in most cases, these 
factors have caused the ratios of the times to peak of the applied 
landing loads to the quarter period of the fundamental mode to approach 
a critical value. 
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Several simplified methods have teen developed for dete rmining the 
inert ia loads in 'wing structures during landing impacts., and studies 
have teen made of the landing-impact inertia loads in simplified 
structures with the use of the principles of these methods. (See refer- 
ence 1 and "bibliography. ) Although experimental investigations have 
been made to determine the magnitudes of inertia loads in actual 
airplane structures (bibliography), little correlation of theory and 
experiment has been made concerning the nature and magnitude of inertia 
loads in airplane -wings during actual landing impacts. 

Another aspect of the problem of wing loads during landing is the 
variation of air load due to changes of attitude and flight path during 
impact . The importance of this change in air load has been the sub Ject 
of some speculation but little investigation. 

In order to evaluate the importance of the various components of 
the load, including dynamic effects and variation in air load, data were 
obtained' during full— scale landing tests of a smaall seaplane to provide 
a comparison of actual wing loads with those predicted by a simplified 
method (reference l). 

The present paper gives a comparison of the theoretical and 
experimental wing loads, in the form of time histories of the wing 
bending moments, and discusses the contributions of each of the components 
of the mcment-( static water— load moment, dynamic water— load moment-, and 
air— load moment) to the total. The static and dynamic components of the 
total moment were calculated and combined according to the procedure of 
reference 1, with the responses of the fundamental mode being calculated 
from the recorded time histories of the applied forces. The air— load 
component was calculated by the procedure of the appendix. 


SIMBOIS 
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lift coefficient / ' \ 

W 

acceleration due to gravity (32.2 ft /sec 2 ) 
lift, pounds 

bending moment in wing, pound— inches 
load factor, multiples of g 
wing surface area, square feet 
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T duration of impact , seconds 

t time, seconds 

t^ time ta peak of applied load, seconds 

t n natural period of fundamental mode, seconds 

u dynamic response factor, ratio of maximum, total water— load wing 

“bending moment to maximum static water— load tending moment 

V Telocity of seaplane, feet per second 

a angle of attack, degrees 

7 flight— path angle, degrees 

A. prefix denoting change 

p density, slugs per cubic foot 

cr response 

cr B static response 

t trim angle, degrees 

co circular frequency, cycles per second 

Subscripts: 
aT average 

c corrected for air load 

h horizontal 

n normal to keel 

o at time of water contact 

p parallel to keel 

r recorded 
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T 

total 

V 

vertical 

max 

maximum 


APPARATUS AND INSTRUMENTATION 


The airplane used in the present . investigation was a small two- 
engine seaplane (fig. 1). Pertinent information about the seaplane is 
given in table I, and additional information may be obtained from 
reference 2. The frequency and shape of the fundamental wing bending 
mode were found from ground vibration tests and are given in table U 
and figure 2. The spanwise weight distribution is also given in 
table H. r ’ "" 

The trim variation was measured with a gyroscopic trim recorder 
mounted in the cabin floor. The airspeed was measured with an NACA 
airspeed recorder, pitot— static— tube' type, . mounted above the cabin. 
Accelerations of the center of gravity were obtained from an NACA 
optical— recording three— component accelerometer mounted securely ln-the 
fuselage near the center of gravity. The time of contact was determined 
from a waten-contact indicator located on the keel at the main step. 

The hull . immersions were determined from pressure gages installed' along 
the bottom of the hull. The wing bending moments were measured by means 
of a strain gage mounted on the wing main span 9 inches from the center 
line of the seaplane (hereinafter referred to as station 9)» 

The estimated accuracies of the experimental data based on 
calibration, instrument, and reading error are as follows: 

Horizontal -velocity, V h> feet per second ±4 

Trim angle, t, degrees . . . . , \. . . . , ... . r . ... . . ±0.25 

Load factor, n, multiples of g ................ ±0’.2 

Initial, wing lift, L 0 ," multiples of— g ...... ±0.05 

Total wing bending moment, Mp, pound.— inches ±0.05 X 10^ 


TEST PROCEDURE' 


The landing-impact' teste, were made in smooth water. During these 
landing tests, airspeed, trim. variation, center— of— gravity accelerations, 
and wing— spar bending momenta were recorded. The landings were made at 
horizontal velocities ranging from 95.4 to 112.0 feet per second, trim 
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angles ranging from 3-00° to 7.83°, and initial flight— path angles 
ranging from 2 . 0 ° to L . b ° . The resulting maximum center— of— gravity 
accelerations normal to the heel line ranged from 1.10g to 1.96g and 
the duration of the impacts varied from O.63 to O.87 second. The times 
to peak of these normal accelerations ranged from 0.10 to O.36 second. 
Values of these parameters and other pertinent information for all the 
tests are presented in table 111. 


THEORY 

Williams 1 Method 


In reference 1 a method was proposed for calculating the dynamic 
effect of an impulsive load applied to an elastic structure. Basically, 
the method follows classical normal-mode vibration theory by considering 
the total response of the elastic structure to a forcing function at any 
instant to be the summation of the responses of all of its normal modes 
at that instant. However , a unique feature, of the method Is that the 
total response of each mode Is separated into .a static and a dynamic 
component,, and the stress due to the sum of the static components of 
the responses of all the modes Is found in one calculation by rigid- 
body analysis. This stress is referred to as the static— load stress. 

The Stress of each mode due to Its dynamic component of response is 
found separately. The total stress Is the sum of the static— load stress 
and the dynamic components of stress for the significant modes. Time 
histories of stress are found in these calculations and thus phase 
relationships of the modes are considered when the stresses for each 
mode are added. 


Air— Load Variation 

Equations are developed In the appendix for determining the effect 
of air— load variation on the wing bending moment during impact. The 
change in bending moment at any instant is expressed in terms of the 
ratio of bending moment to lift at time of contact and the wing lift at 
that instant. In developing the equations, the air load Is assumed to 
change Instantaneously with change In angle of attack and the rate of 
change in air load Is assumed to be slow enough to neglect structural 
dynamic response. Also, the ratio of bending moment to lift Is assumed 
to be constant throughout the Impact. 
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CALCULATIONS AND RESULTS 


The variations of the ylng angle of attack and velocity during 
impact necessary for computing the changes in. tending moment by 
equation (4) of the appendix were determined for each- impact from the 
recorded data in the following manner. The accelerations normal and 
parallel to the keel line, obtained from the three— component center-of— 
gravity accelerometer* were plotted (fig. 3). After the trim-angle 
variation (fig. 4(a)) was taken into account, these accelerations were 
resolved into vertical, and horizontal components. Integration of the 
time histories of these accelerations over the duration of the Impact 
produced time histories of the changes in vertical an<J. horizontal 
velocities. Since' the vertical velocity at the time- of contact was 
not accurately known, the initial velocity was determined so that 
integration of the time-history curve over the duration of the impact 
resulted in a final vertical displacement of zero. (The duration of 
the impact is defined as the interval between the time of ^contact and 
the time when the center of gravity again reached Its initial height 
above the mean water line. The Instant of contact was found from a 
water— contact indicator on the step and the time history of the center- 
of— gravity displacement was determined from the times of limners ion and . 
emersion of the hull pressure gages, the fixed location of the center 
of gravity relative to the step, and triut-angle time history. ) 

Integration of the time— history curve of the corrected vertical velocity 
from time of contact to the time of zero vertical velocity determined 
the maximum displacement of the center, of gravity. The maximum 
displacements determined In this manner for all the impacts agreed 
within experimental error with the ma-H mum displacements calculated from 
the hull pressure gages. With the use of the corrected vertical— velocity 
and horizontal— velocity time histories, time histories of the flight— path 
angle y and the resultant velocity were computed. From the time 
histories of trim angle t (fig. 4(a)) and flight— path angle y 
(fig. 4(b)), the time history of the angle of attack a was computed 
(fig. 4(c)). 

With the . use of the time histories of angle of attack and resultant 
velocity, the changes In bending moment In the wing at station 9 due to 
the changes in air load were determined for each impact by use of the 
equations in the appendix and are presented in parts (a) of -figures 5 
to 10 . 

The procedure of reference 1 was used to compute the bending 
moments because it provides a convenient means of applying the principles 
essential to a dynamic— loads analysis which results In time histories 
of the wing bending moments. The forcing function for each' impact was 
determined from the normal acceleration measured in the hull by an 
accelerometer located near the seaplane center of gravity. Because of 
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the manner of* connection "between the wing and hull (fig. l), the 
measured acceleration was not appreciably affected "by the oscillations 
of the ~w ing . By including the effect of the varying wing lift on the 
center— of— grarity acceleration,, the acceleration normal to the keel 
line due to the hydrodynamic force only was determined. Prom this 
acceleration and the mass of the seaplane, the hydrodynamic forcing 
function was calculated. The dynamic responses of the significant 
inodes to this forcing function were computed "by a recurrence method 
developed at the Langley Structures Research Division in which the 
actual forcing function is used without approximation. Only the 
dynamic effects of the fundamental’ "bending mode were included in the 
final results "because calculations showed the dynamic effects of the 
second symmetrical bending mode to be negligible. This observation 
was borne out by the absence of higher mode effects on the strain— gage 
records. The calculated time histories of the Btatic water-load and 
dynamic water— load components of wing bending moment at station 9 are 
presented in parts (a) of figures 5 "bo 10. The s panwise bending- 
moment distribution for the fundamental-mode 1 g inertia loading 
calculated as set forth in reference 1, a 1 g static water loading, and 
a level— flight loading are plotted in figure 11. The values of the 
bending moment at station 9 used in the application of the method of 
reference 1 were obtained from this plot. 

These three components of wing-bending-moment changes obtained in 
this manner for station 9 were combined and added to the wing bending 
moment existing' at the Instant of contact. This total theoretical wing 
bending moment is presented in parts (b) of figures 5 "to 10 together 
with the wing-bending-moment variation measured by the strain gage at 
station 9. 

In order to demonstrate the accuracy gained by using the actual 
forcing function, responses were computed for two typical Impacts with 
the use of apparently good approximations to the recorded forcing 
functions. Parts (a) of figures 12 and 13 show the actual forcing 
functions for these Impacts together with the approximations. The 
responses to the approximations are presented in parts (b) of 
figures 12 and 13 together with the responses obtained from the actual 
forc in g functions with the use of the recurrence method. 


DISCUSSION" 


Comparisons between the total theoretical and experimental wing 
bending moments are presented In parts (b) of figures 5 to 10. The 
comparisons are made only for wing station 9 because the bending moments 
in the outer wing section were so small as to be of the same order as 
the estimated error. Only the dynamic effects of the fundamental 
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■bending mode vere included in the final results "because calculations 
shoved the dynamic effect's of the second symmetrical "bending mode to he 
negligible. This negligibility vas borne out by the absence of- higher- 
mode effects on ’the strain— gage records. The comparisons shov the 
predicted values, to be in good agreement with the experimental values. 

As can be seen from the figures, the phase .3^elation§hips betveen the 
theoretical and experimental values are consistently good, and' the 
maximum changes from initial conditions shov a range of error of 5 
to 28 percent based on the experimental values of the maximum changes 
in ving moment. These results . indicate that vhen the three components 
of moment are included in the theory, good agreement is obtained, 

A comparison of the level— flight bending moment and the change in 
bending moment . due to the heaviest impact-of the tests shovs that- the 
maximum change ~ in bending moment- accompanying a dovnvard motion of the 
ving vas approximately 50 percent of the level— flight moment (figs. 11 
and 5 (h))- This maximum bending-moment change vas produced by a 
2 . 03 g impact, . When any differences in the response factor and any 
change in air— load bending moment- are neglected, landings of over 4 g 
vould be required for the dovnvard motion of the ving merely to start 
stress reversal in the ving. Similarly, even if the maximum dynamic 
vater— load bending moment (fig. 8(a)) caused by a 1.9^g impact vere 
to be tvi.ce as large for a i+g impact, vere to exist after the vater' 
load vas removed, and vere to be superposed on a level— flight moment— 
the moment produced by the upvard motion of the ving vould still be 
less thantvice the level— flight moment. Therefore, the- change in 
bending moment hue to a landing impact is unimportant in this seaplane 
insofar as this change vill not produce critical stresses at the ving 
root. This unimportance may be largely attributed to the fact- that the 
fundamental-mode 1 g' inertia loading is relatively small as compared 
vith a level— flight - loading (fig. 11). 

The effects of structural flexibility on the computed bending 
moments can be seen by comparing the static and dynamic components of 
the vater— load bending moment - in parts (a) of figures 5 to 10 and by 
observing the dynamic response factors u in table XII. The dynamic 
overstresa attributable to structural flexibility is the dynamic 
component of stress -In parts (a) of figures 5 to 10 and is represented 
in the response factor u by the amount thaf this factor differs from 
unity. Since re sponse_f actors as high as 1.97 &re obtained, the dynamic 
overstress sometimes contributes an increment of stress almost as large 
as the static vater^-load stress. This observation is in agreement vith 
the results of other investigators (bibliography) and shovs the necessity 
of using dynamic analyses in landing— load investigations. 

The change in ving bending moment due to change _ in air load on the 
ving is a function of the changes in velocity and angle of attack. (See 
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the appendix.) In these tests the changes in Telocity during the 
impacts were small (table HI) and the. changes in airload bending 
moment vere therefore almost entirely attributable to the changes in 
angle of attack. ' Since the angle of attack is a function of the trim 
and flight— path angles, large changes in angle of attack will occur 
when the trim and flight— path angles change to a great degree. For 
the relatively small changes in trim angle and flight— path angle which 
occurred in these tests (the ma xim um values being 9*31° and 7*0 ° s 
respectively), air-load changes as large as 0 . 2 g were computed. For 
this airplane, these changes in air— load bending moment were of about 
the same magnitude as the b ending-moment changes due to the static 
component of loading (parts (a) of figs. 5 to 10) and inclusion of the 
effects of air— load changes in the calculations was therefore necessary. 
For other airplanes with structural and mass characteristics conducive 
to larger inertia— load moments, the bending-moment changes accompanying 
a 0.2g change in air load would be small relative to the changes in 
bending moment caused by inertia loads. However, for more severe 
changes in flight— path angle, which should be considered in a design 
analysis, the effects of the change in air load on the bending moment 
in the wing may still be large enough to warrant consideration in a 
design analysis. Further investigation is necessary to determine the 
importance of this air— load variation in design— strength calculations. 

In most landing tests the applied forcing functions are not 
easily obtained from the center— of— gravity accelerations because of the 
superposed accelerations caused by structural oscillations. But 
because the fuselage of this airplane is connected to the wing by struts 
located near the nodal point of the fundamental mode which represented 
the greatest portion of the wing bending and by a nonstructural fairing 
which neither transmitted nor interfered with the wing oscillations, 
the accelerometer in the fuselage was not appreciably affected by the N 
wing oscillations in these tests. 

The calculation of the dynamic response of each of the normal modes 
of the seaplane involves solving for the response of an equivalent 
simple mass oscillator to the given forcing function. In this 
calculation it is common practice to approximate the forcing function 
in order to simplify the computation. In order to indicate the 'effect 
of such approximations on the accuracy of the response, figures 12 
and 13 have been prepared. ■ In these figures are shown two typical 
forcing functions, various simple approximations to each, and the 
responses computed from the two forcing functions and their approxi- 
mations. By comparing the calculated responses of the approximations 
to those of the recorded forcing functions in figures 12(b) and 13(b), 
it is seen that errors as large as 2 Q percent in the total response 
may be introduced into the calculations by use of apparently good 
approximations. A recurrence method developed by the Langley Structures 
Research Division was used to compute the total responses to the recorded 
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forcing functions used in proparing figures 12 and 13 . In this 
procedure the motion of each of the significant normal modes of the 
airplane is computed from the equation of motion of„an equivalent 
simple mass oscillator subjected to the recorded forcing function. 
Figures lh(a) and ld(b) have been prepared to show. the accuracy with 
which the total response of the equivalent simple mass oscillator can 
be obtained by this procedure. Comparison Is made of the ratios of 
the total response or to the maximum static response of the normal 
mode o'smax for siri 5 ) - 1 - e sinusoidal and triangular forcing functions 

calculated by both this method and Duhamel 1 s integral. It can be seen 
that excellent^ agreement is obtained. 


CONCLUSIONS 


Experimental wing bending moments obtained from a landing— impact 

Investigation with a small full-scale seaplane were ..'.compared with 
analytical results. The effects of the variation in air load during 
the impacts were included in the analytical procedure. The responses 
of the fundamental mode were calculated from the recorded time histories 
of the applied hydrodynamic forces. For the seaplane tested and the 
conditions of the Impacts encountered, the following conclusions may be 
drawn: 

1. Good agreement between measured and computed wing bending 
moments was obtained when the three components of the wing moment 
(static water— load moment , dynamic water-J.oad moment, and air— load 
moment) were Included in the calculations. 

2. The effects of structural flexibility on the wing bending 
moments, represented by the dynamic overstress and understress, were 
large, the moment- -due to the dynamic component of the total response 
being as much as 97 percent - of that caused by the static water-load 
component. 

3 . Although the changes in seaplane attitude during the landing 
impacts were small, the variation in the air-load component of the 
total moment was of about the same magnitude as the. .static water— load 
component. Although this comparison of changes Is not representative 
of the relative importance of the air— load variation In seaplanes with 
structures conducive to large Inertia— load components. It indicates 
the probable significance of the effects of air— load variation since 
large changes in seaplane attitude must - also be considered. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va., December 21,.. 19^9 
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APPENDIX 


DETEHMINATION OF THE EFFECT OF CHANGE IN AIR 
LOAD ON WING BENDING MOMENT 


In order to determine the effect of change In air load on the 
"bending moment in the wing during impact, a procedure -was developed to 
determine the moment at any time t in terms of the ratio of "bending 
moment to lift at time of contact and the wing lift at time t. The 
ratio of "bending moment to lift was assumed to "be constant, or 


M 

L 



Therefore , 



Io 


or 


M, 


M = — i pS V 2 ^ 


L o 2 


The ratio Mq/Lq was obtained from initial conditions. The change in 
moment may be expressed as 


M o 1 2 

M° " “ = “ g P^o C] 


- - - poT^C-, 
L 2 ^ 3 
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or 


AM 



V o 2 c Lo/ 


( 1 ) 


The numerical value of may Ido calculated from the lift equation 

and Initial conditions. Therefore, 



h> 


or 



where n Is the acceleration In multiples of g acting on the seaplane 
at contact and W ' is the seaplane weight - in pounds. 

From the curve of plotted against a (no flaps) supplied "by 

the manufacturer, the initial conditions were found to lie in the 
straight portion of the curve for all impacts considered. The equation 
for . Cl may therefore "be expressed as 7. 


Ct — t 4- m co 


( 2 ) 


where K is the Cl intercept - at a = 0° and m is the slope of the 

curve. The value of K was determined from the initial conditions of 
the impacts. Since the slope of the curve of Cl plotted against ex 

is constant regardless of flap setting, an average value of £ was 
"believed to he .usable, provided the error in calculating the lift with 
the use of this average value would not he greater than that estimated 
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for the recording instruments. Accordingly ^ such, a procedure was carried 
out and the maximum, error was 3*6 percent. Since the slope of the lift 
curve is known to he 0 . 014 - 98 * equation ( 2 ) may he written 


C L * ^av + °* oi *-98a. 


(3) 


Suhatitution of equation ( 3 ) in equation (l) for Cp results in 


Mq 1 o 

*** - £ ! ^ % 


1 - 


T 2 (Ea T + 0.0498a) 


T o\ 


or 


m = b 


1 


y 2 (k^ t + 0. 0^-980,) 

c. 


w 


where 


_ Mo 1 , 2.. 


o 


C = 

o 


Equation (4) eatresses the change in bending moment £> M due to 
change in air load in terms of initial conditions and measured 
variables . 
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TABLE X 

GENERAL UXFORMATIOF ABOUT SEAPLANE USED IN 
LANDING— IMPACT INVESTIGATION 

Approximate flying weight Turing tests, It ...... 19,200 

Stalling speed (flaps down), fps 9^ 

Wing area, sq. ft. ...<,. j8o .6 

First natural frequency, cps . . k.j6 

Second natural frequency, cps 13 .0 
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TABLE II 

WEIGHT DISTRIBLTION AND MODE-SHAPE FACTOR FOR FUNDAMENTAL 
BENDING (L.T6 cps) MODE OF WING SEMISPAN 









TABUS m 

TT.TWT AHD IMPACT YAKEABLE3 
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1. NACA optical-recording three - 

component accelerometer 

2. Gyroscopic trim recorder 


3. Pressure gages and water contacts 
Water -contact indicator 


Figure 1.- Seaplane and instrumentation. 


Mode-shape factor 


i 



Distance from seaplane center line, in. 

Figure 2.- Experimental fundamental bending-mode shape of wing semispan. 

Frequency, k.’jG cycles per second. 
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Acceleration, multiples of g 
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Figure 3 .- Typical accelerometer records a&d fairings. 
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(a) Trim-angle variation during Impact. 
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(b) Flight -path-angle variation during impact. 





Time, sec 

(c) A ngle-of -attack variation during impact. 
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Wing "bending moment, lb-in. 
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Figure 10.- Wing -tending -moment time histories during impact, run 
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Figure 11.- Level-flight 
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Level -flight Lending moment 

1 g static mb. ter -load.. "bending moment 

Fundamental-mode 
1 g inertia -load Lending moment 


Forcing function, n, multiples of g 
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(a) Recorded forcing function and approximations. 




Figure 12.- Recorded forcing function for run 4, approximations, and calculated responses. 
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0 .1 .2 .3 . 1 * .5 .6 .7 

Time, sec 

(b) Calculated responses to forcing function and approximations. 


Figure 12.- Concluded. 
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(a) Recorded forcing function and approximations. 

Figure 13 *- Recorded forcing function for run 5 * approximations, and calculated responses. 
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0 .1 »2 *3 »4 >5 »6 »7 

Time, sec 

(b) Calculated responses to forcing function and approximations. 


Figure 13.- Concluded. 
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(a) Responses to sinusoidal forcing function. 


Figure 1^.- Comparison of total responses. 
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(b) Responses to triangular forcing function. 
Figure l4.- Concluded. 
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